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c§ ■ Abstract 

The description of an impurity atom in a Bose-Einstein condensate can be cast in the form of 
^ , FrohHch's polaron Hamiltonian, where the Bogohubov excitations play the role of the phonons. 

' An expression for the corresponding polaronic coupling strength is derived, relating the coupling 

qh' strength to the scattering lengths, the trap size and the number of Bose condensed atoms. This 

allows to identify several approaches to reach the strong-coupling limit for the quantum gas po- 
larons, whereas this limit was hitherto experimentally inaccessible in solids. We apply Feynman's 
path-integral method to calculate for all coupling strengths the polaronic shift in the free energy 
' and the increase in the effective mass. The effect of temperature on these quantities is included in 

Q ' the description. We find similarities to the acoustic polaron results and indications of a transition 

, between free polarons and self-trapped polarons. The prospects, based on the current theory, of 

investigating the polaron physics with ultracold gases are discussed for lithium atoms in a sodium 
^ ! condensate. 
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I. INTRODUCTION 



Quantum gases are used as an excellent testbed for many-body theory, and are particu- 
larly useful to investigate strong-coupling^ regimes or strongly correlated regimes that have 
remained out of reach in the solid state[l]. The current work focuses on the physics of an 
impurity in a Bose-Einstein condensate (BEC), as this system opens a new and promising 
avenue to investigate long-standing problems in polaron theory [3]. In 1952 Frohlich, inspired 
by Landau's concept of the polaron, derived a Hamiltonian that describes a charge carrier 
(electron, hole interacting with its self-induced polarization, in an ionic crystal or a polar 
semiconductorjl, 0]. The Frohlich polaron Hamiltonian has resisted exact analytical diago- 
nalization since 1952 and became a laboratory to test methods of quantum field theoryj^. 
Tomonaga's canonical transformation was applied to study the weak coupling regime. Bo- 
goliubov tackled the polaron strong coupling limit with one of his canonical transformations. 
Feynman used his path integral formalism (and an ad-hoc variational principle) to develop 
a superior all coupling approximation for the polaron [5i]. These studies were addressing the 
self-energy, the effective mass and the mobility of Frohlich polarons. Studies of the optical 
absorption of a Frohlich polaron were initiated by Evrard et al. for strong coupling and 
by Gurevich et al. \j\ for the weak coupling limit. Devreese et al. calculated the op- 
tical absorption of the Frohlich polaron at all coupling using Feynman path integrals and 
the response formalism introduced in Ref. [9|. The results of Ref. jsf reveal the internal 
excitation structure of the Frohlich polaron that is absent in Ref. [9|. 

However the largest Frohlich polaron coupling constant established in any solid by ex- 
periment is not large enough to reveal the rich internal excitation structure of the optical 
absorption predicted in Ref. jl]. The possibility to tune the coupling strength in the BEC- 
impurity system presents a marvellous challenge and opportunity to experimentally reveal 
the internal excitation structure and its related resonances and scattering states "contained" 
in the Frohlich field-theoretical polaron Hamiltonian. Quantum gases could therefore reveal 
important and subtle characteristics (eigenstates, resonances,...) of a Hamiltonian devised 
for a solid, however, that cannot be realized in a solid. The observation of the spectra of 
the BEC-impurity Hamiltonian at large coupling would be all the more interesting in view 
of recent studies of the spectrum of the Frohlich polaron Hamiltonian with Diagrammatic 



Quantum Monte-Carlo numerical techniques [10| and with a new strong coupling model [11 



that explore the validity of the Franck-Condon principle with increasing coupling strength. 
Similarly, the transition between quasi-free polarons and self-trapped polarons predicted for 
acoustic polarons[li] is expected to occur in a coupling regime hitherto inaccessible in the 
solid state. Finally, a better understanding of the intermediate and strong-coupling regimes 
is needed to further elucidate the role of polarons and bipolarons in unconventional pairing 
mechanisms for high-temperature superconductivity [l^ . 

The Hamiltonian of an impurity in a Bose-Einstein condensate can be mapped onto the 



polaron Hamiltonian when the Bogoliubov approximation is valid[l5|, [IGJ. The polaronic 



effects comes about through the coupling of the impurity with the Bogoliubov excitations, 
as shown in Sec. II. We apply the path-integral formalism, valid at all values of the coupling 
strength, to describe the polaronic effect in a condensate in Sec. HI. As a concrete example, 
we use the parameters for a lithium impurity in a sodium condensate, and derive an expres- 
sion linking the experimental parameters to a dimensionless coupling constant analogous to 
the Frohlich coupling constant for electrons in polar crystals. The polaronic energy shift, 
the effective mass increase and the polaron radius are studied for all values of this coupling 
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strength and at all temperatures where the Bogoliubov approximation is valid, in Sec. IV. 



II. THE POLARON HAMILTONIAN FOR AN IMPURITY IN A CONDENSATE 

The Hamiltonian of a single atom, in the presence of a Bose gas is given by: 



f) 1 

H = ^ + £k4«k + 2 5Z ^BB(q)4'-q«k+q«kak' + ^/i?(q)p/(q)4'-q«k' (i) 

k k,k',q k,q 

The first term represents the kinetic energy of the "impurity" atom with mass rrij. The 
operators a|^,ak create and annihilate a boson with mass ms, wave number k, and energy 
£k = {hk)"^ / {210113) — n where /i is the chemical potential. These bosons interact, and Veslq) 
is the Fourier transform of the boson-boson interaction potential. The interaction between 
the bosonic atoms and the impurity atom is described by a potential Vib{^) coupling the 
boson density to the impurity density pi{q) which can be expressed as a function of the 
impurity position operator f as 

p/(q) = j dh e-''i-'"5(r-f). (2) 

Bose-Einstein condensation is realized when the single-particle density matrix has an 
eigenvalue Nq ^ 1 comparable to the total number of bosons [l^]. For a homogeneous 



condensate, this is expressed by the Bogoliubov shift |20l] which transforms the Hamiltonian 
(P) into 



H = Egp + NoVjBiSi) + ^ + J^^k&i&k + 5^ J^V^/B(k)p/(k) [Sk + hl^ (3) 

kT^O k^O * ^ 

Here, the operators h\., 6k create resp. annihilate Bogoliubov excitations with wave number 
k and dispersion 

£^k = A/ek[ek + 2iVoVBB(k)] (4) 
where ^k = (^fc)^/(2mB). The first term in (j3]) represents the Gross-Pitaevskii energy [2l[ 

2 ■ 2 

k^O 

and the second term is the interaction shift due to the impurity. For both the boson-boson 
interaction and the boson-impurity interaction we will assume contact pseudopotentials: 
VBs(r — r') = (yfBB5(r — r') and ^/^(r — r') = gjBSij — r'). The interaction strengths Qbb and 
giB are related to the boson-boson scattering length a^s and the boson-impurity scattering 
length aiBi respectively, through the Lippmann-Schwinger equation. For the boson-boson 
interaction, the first order result qbb = ^'nf^aBB I'^b will suffice. For the boson-impurity 
interaction, the Lippmann-Schwinger equation needs to be treated correctly up to second 
order to obtain valid results for the polaron problem, since as we shall see, (^f^ will appear 
in the expressions. 



AJ2 1 

Egp = NoBo + -fVBBiO) + -YNoVbb{^) (5) 
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The resulting Hamiltonian ([3]) maps onto the Frohhch polaron Hamiltonian 
with 



k^^O 



kT^O 



(6) 



1/4 



giB 



(7) 

(8) 



In these expressions ^ = 1/ vSvraBsno is the heahng length of the Bose condensate with uq = 
Nq/V the condensate density (in the present calculatations for the homogeneous gas we work 
with unit volume) and c = h/ (v^m^^) is the speed of sound in the condensate. The operator 
structure of the BEC-impurity Hamiltonian and that of the Frohlich polaron Hamiltonian 
is identical. Depending on the analytical form of the scalar functions multiplying the field 
operators, the Frohlich polaron Hamiltonian, originally devised to describe the electron/hole 



- longitudinal optical phonon interaction, can depict as well the acousto-polaron[l^, Il3| . 
the piezo-polaron[l3], the ripplo-polaron [l8|, etc... However, a formalism/method that 
diagonalizes the Frohlich Hamiltonian at all coupling, does so for all those different types of 
polarons, including the BEC-impurity polaron. 

This mapping of the BEC-impurity problem onto the Frohlich Hamiltonian, and the re- 
sulting possibility of polaronic self-trapping of an impurity in a condensate, has been the 



subject of recent theoretical investigations at strong coupling [15|, [IGj, and at low coupling 
with the Lee-Low-Pines scheme 22|]. For small polarons in an optical lattice, a perturba- 
tive study exists [2^. As is also known from the study of the polaronic problem for slow 
electrons in a polar crystal, these approximations give qualitatively different results for dif- 
ferent regimes of coupling strength. These different characteristics between weak-coupling 
and strong-coupling behavior are even more pronounced in the case of an impurity in a con- 
densate. In this work, we present a calculation valid for all coupling strengths, based on the 
Feynman- Jensen variational scheme first successfully applied by Feynman ^] to study slow 
electrons interacting with optical phonons in polar crystals. As discussed below, the present 
system will have more similarities to acoustic polarons than to the Frohlich polarons. 



III. JENSEN-FEYNMAN FREE ENERGY 
A. General treatment 

The calculation of the free energy F of the polaron is based on the Feynman- Jensen 



variational inequality [2J, |25 



F^F, + ^{S-So),^ (9) 



where S is the action functional of the real system described by the Hamiltonian ([6]), and So 
is the action functional of a variational model system with free energy Fq, and (3 = l/(/csT) 
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with T the temperature. This variational principle is an extension of the usual Gibbs- 
Bogoliubov variational principle F ^ -Fo + {h — HqJ where H is the Hamiltonian of the 

system under study and Hq is a model Hamiltonian, and the brackets indicate statistical 



averaging. The inequality [9] can be derived 2J, |25|] from the path-integral expression for the 
partition sum: 

Vr e-{'5[r(*)]-5o[r(t)]}/ftg-5o[r(t)]/ft 

= e"^^" /g-(5-5o)/ft\ 

\ / So 

^ exp [-/?Fo -{S- So)s, /h] (10) 

where the inequality follows from the convexity of the exponential, and Vr represents the 
path-integral measure. The peculiarity of this variational principle, in comparison to the 
usual Rayleigh-Ritz variational approach, is that a variational action functional is used 
rather than a variational wave function. One important strength of the Feynman- Jensen 



principle is that it can be used at nonzero temperatures [2J]. Another strength, in comparison 
with the Gibbs-Bogoliubov variational principle, is that it can treat systems that are hard 
or impossible to describe in the standard Hamiltonian formalism, such as a collection of 
particles that interact through a retarded potential, such as Lienard-Wiechert potentials 
in electrodynamics, or the phonon-induced electron-electron potential in polaron theory. 
In particular, the polaron Hamiltonian ([6]) after the elimination of the phonon degrees of 
freedom leads to an action functional S containing retardation effects: 

ft/3 2 

S = j '^f^r)dr - Y}^^ jdrjda ^(k, \r - a\)e^^-^^^^^-<'^\ (11) 



with Q (k, u) the Bogoliubov excitation Green's function, given by 

r(\. .,\ cosh [u)^{u - hl3/2)] 

^('''")= sinh(n/5c.,/2) ' (^2) 

The action functional (fTTj) is obtained by integrating out analytically the oscillator degrees 
of freedom corresponding to the Bogoliubov excitations, a technique introduced in the con- 
text of phonons in Ref. [5|. This gives rise to a retarded interaction (proportional to 
gik-[r(r)-r((T)]-j^ mediated by the Bogoliubov-excitation Green's function, and with coupling 
strength \V]^\^ /{2h). 

The system under study here is modelled by a variational trial system with free energy 
Fq described by the action functional 

ft/3 ft/3 ft/3 

fmi.^,^^ MW^ f , r cosh[W{\r-a\-h(3/2)] ^ ^ ^ ,,^2 
5„ = y -rHr)dr + —J ir j ia^-^L_^^ |r(r) - r(a)f . (13) 



The model system corresponds to a mass m/ that is coupled by a spring with spring con- 
stant MW'^ to a "Bogoliubov" mass M. Both M and W are variational parameters. The 
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expectation value in the inequality is evaluated in the model system. This model system is 
chosen partly because the expectation values can be calculated analytically for it. For more 
details on the formalism, we refer to Ref. [sl, [13, 26| : 



F ^ 



/5 
3 



In 



2 sinh 



-In 



2 sinh 



(^) 



M 



2(3mi + M 

hf3 



h 



du 1 



u 

h(3 



2 J 

g{k,u)MM,n{Ku) 



(14) 



Here Q = Wy^l + M/rrij, and the memory function is given by 



MM,n{Ku) = exp 



2{mi + M) 



u 



u — 



M cosh [Vt{hl3/2)] - cosh [Vt{hl3/2 - u)] 



h(3 mi 



VLsi-nh{l3nn/2) 



15) 

Rather than choosing VF, M as variational parameters, we can choose f2, M as variational 
parameters. 



B. BEC-impurity 



Care should be taken when we substitute the interaction amplitude ([HD into these expres- 
sions. Since |Vk|^ oc gj^, we need to solve the Lippmann-Schwinger equation up to second 
order to obtain the link between qib and ajB correctly; we find 



rrir 



(16) 



where is the relative mass (m^ ^ = + mj^). From this, we find that the term 
NqVib{0) = Nog IB in ([3]) must also contribute a renormalization factor since 

NogiB No 



rrir 



(17) 



Substituting ([8]) and into ([HD then yields 



F<-Un 

3 M 



2 sinh I — 



pn 



In 



2/51 + M 



^^^^ — coth ( ^— 



2 sinh 
1 



pn 



2vrTM 



(18) 



+ a 



mB + 1 
rfiB 



dk 
An 



2m 



B 



/9/2 



l+mBj V2TP 



Q{k,u)M{k,u)du 
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In this expression, 

« = O/i?/(«BB0 (19) 

is the polaronic couphng strength. We use polaronic units, so that H = mj = ^ = 1. That 
means the free energy is in units h"^ /{mi^'^), so that (3 = h'^/{mj^'^kBT) and rfiB = rriB/mi. 
The integration variables k, u and the variational parameters M, VL in f[T5]) are dimensionless. 
In these units, we find after substitution of ([H]) and ([7]), that the Green's function is given 

by 

^ cosh [kVWT2{2u - (3)1 {AmB)] 
~ sinh [/5fcv/FT2/ {AmB)] 
and the memory function takes the form 



J^{k, u) = exp 



fu{f] -u) ^ ^cosh [Q{f3/2)] - cosh [Q{(3/2 - u)] 



2(1 + M) V f3 nsmh{(3n/2) 



(20) 



The first term of the integrand of the fc-integral in (fT8|) is due to a renormalization factor 
which arises when we relate qjb to the boson-impurity scattering length. This term is 
independent of the variational parameters M and Q and does not infiuence the optimization 
of the variational parameters; it is however necessary for the /c-integrals to converge. In the 
case of the acoustic polaron, convergence in the variational minimization of (fTSl) is sped up 



by introducing a cut-off Kc to the wave-number integral [12|, Il3|, for which there is a natural 
choice, namely the edge of the Brillouin zone. For atomic gases, as we show in Sec. IV, the 
results do not strongly depend on the value cut-off as long as it is on the order of the (inverse) 
Van der Waals radius which we use as a natural scale for Kc in the present case. Note that, 
although these formulae are derived for any temperature (any (3), the mapping onto the 
polaron Hamiltonian is only valid for temperatures below the critical temperature, and low 
enough so that the number of Bogoliubov excitations can be neglected with respect to the 
number of atoms in the condensate, so that ([3]) gives a good description of the quantum gas. 
In the discussion, we will look at realistic experimental values for (3 and other parameters. 



C. Weak and strong coupling limits 



From our central result flTSl) . we can retrieve both the weak-coupling and the strong- 
coupling results by a judicious choice of variational parameters, or rather variational trial 
actions. The weak-coupling limit is obtained by taking the limit M — >^ in the Jensen- 
Feynman treatment, which corresponds to a the trial action 



hf3 



^weak 



mi .2 



r (T)dT 



(21) 



of a free particle, as pointed out in the original work by Feynman [5j. In the limit of zero 
temperature, this corresponds exactly to the second-order perturbation result E^^"'^{a); we 
find 



E 



weak 



a- 



rriB) 



rriB 



2tx 



dq 



l+rfiB Y g2 -^2mBq + \/ + 2 



(22) 
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Note that the variational energy becomes independent of Q in the hmit M ^ 0. The energy 
grows hnearly with a, and the proportionahty constant depends on : 

{a /{V27c) for rriB oo 
4V2a/(37r) for 1 (23) 

a /(v^ms) for rriB 

The strong-couphng hmit is derived by taking the hmit M ^ oo in the Jensen- Feynman 
treatment. In this hmit, the variational particle with mass M — > cxd is fixed in place, and 
we get a particle attached by a spring to a fixed point rather than to a mobile variational 
mass. The resulting trial action for strong coupling is[5| 

(- ■ 

2 ^ ' 2 



gstron, ^ I ^ "^T^t) ) dr. 



relating the strong-coupling limit to the variational approach of Landau and Pekar[27| for 
the polaron problem, based on a Gaussian variational wave function. The variational energy 
becomes 



^strong < _^ + 



2 °° 

3 /mB + l\ [ dk 



a 



4 p^oo \ niB J J Att 







(3/2 

2m R \ k^ 



l + rfiBj V2T¥ 



Q{k, u)M.{k, u)du 







(24) 

Care must be taken with the limit (3 ^ oo : since the convergence of Q{k, u) to its /? — > cxo 
limit is not uniform the limit can not be interchanged with the integration. 



IV. RESULTS AND DISCUSSION 
A. Experimental realization 

Before we turn to the results, it is useful to find typical experimental values of the system 
parameters. Consider a sodium condensate of A^o = 10^ atoms in a harmonic trap with 
axial and radial trapping frequencies c^ax = 27r • 50 Hz and tUrad = 27r ■ 150 Hz, respectively. 
As scattering length for sodium in the \F^Mp) = |1,1) hyperfine state, we use aBB = 2.8 
nm 28|. The condensate is suitably described by the Thomas- Fermi (TF) approximation. In 
our calculations we have assumed spatial homogeneity. The density appears in the theoretical 
results through the healing length, which is the relevant length-scale for the polaronic effect, 
and appears in expression (fT9l) for the coupling constant. As long as the polaronic effects 
are restricted to a range considerably smaller than the TF radius (9.5 /xm), we can use the 
local density at the center of the trap, hq = 7 x 10^^ cm"^, to estimate the appropriate 
healing length. In the current example, it is ^ = 450 nm at the center of the trap, indeed 
smaller than the TF radius. Using ^ = l/^/SirnoaBB and the Thomas- Fermi expression for 
no in the expression fll9p . we find the useful relation for the coupling constant 



1/5 

V dffo J ctBBdno 



a = 15'/'(^^) (25) 
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expressing the polaronic coupling constant as a function of the number of condensate atoms 
A^o, the oscillator length ano = a//i/ (mujHo), associated with the geometrical average trap 
frequency ujho = ^/^ax^rad^ the s-wave scattering lengths asB and ajB- 
For the impurity, we consider a ^Li atom, so that 

niB = mB/rrij = 3.8 (26) 

Using ajB = 0.8 nm as the ^Li-^^Na scattering length 29|, a is of the order of 10^'^. For 



these values, the weak-coupling expression fl22|) is suitable. To increase a, several strategies 
are possible. One possibility is to increase ajB/o^BB using Feshbach resonances to either 
increase ajB or decrease abb. One could also increase aiB/ano by tightening the trap. 
The remaining factor in expression fl25|) varies slowly as a function of the TF parameter 
NoaBB/aHO- 

We look in particular at two strategies for increasing a in the case of lithium impurities 
in a sodium condensate. On the one hand, one can minimize qbb using the zero-crossing 
near the sodium Feshbach resonance at 907 GjSOj. The tunability of a^e is limited by the 
magnetic field stability; we expect a factor of 10 increase in a to be realistic. On the other 
hand, one can increase qib using a sodium-lithium Feshbach resonance, for which the prime 
candidates are the already observed resonance at 796 G [31] or the predicted resonance at 
1186 G [29!]. In this strategy, the limiting factor to the tunability of ajB is the trapping 
lifetime of the lithium atoms. The main loss process is three-body recombination, involving 
one lithium and two sodium atoms, which is strongly enhanced near the lithium-sodium 
resonance as the loss scales as a^^. We expect that an increase in a (oc a]^) up to a factor 
1000 are still attainable. This leads to a maximal a of order unity, which is in the crossover to 
the strong-coupling regime. In addition, to avoid phase separation or collapse the trapping 
potential and the number of lithium atoms have to be carefully chosen. 

The energy scale ^^/(m/^^) corresponds to a frequency of 27r ■ 8.1 kHz or a temperature 
of ca. 390 nK. With this choice, 

B = (27) 
^ T (nK) ^^^^ 

The (ideal gas) critical temperature for the sodium condensate is 220 nK, 01 (3 = 1.8. As 
the lowest temperature we estimate 0.5 where = n/kB = 53 nK. This corresponds to 
(3 = 15, so that the experimentally relevant window for this parameter is /? G [1.8, 15]. 

The Bogoliubov dispersion at low k is similar to that of acoustic phonons rather than to 
that of optical phonons. For the acoustic polaron, it is known that the variational parameters 
-especially the polaron mass- depend on the value of a cut-off in the integralpJl]. For 
phonons in solids, this cut-off is related to the edge of the Brillouin zone. In the case of dilute 
quantum gases, a natural cut-off scale arises from the range ro of the interatomic potential: 
on scales smaller than this the interaction amplitude cannot be represented any more by 
expression (IHl) and should be suppressed. The range of the interatomic potential for sodium 
is estimated through the Van der Waals radius tq = 2.4 nm, related by tq = ^(mCg/Zi^)^/^ 



to the Van der Waals coefficient Cq = 1556 |32|. In units of ^ ^, this places the cut-off at 



Kc ~ 200 for the parameters listed in this section. 



B. Free energy and critical coupling strength 

Figure [1] shows the results for the free energy as a function of a, using mB/mj = 3.8 
and a wave number cut-off at Kc = 200, for different values of /3, ranging from /3 = 100 to 
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FIG. 1: (Color online) The polaronic contribution (jlSp to the free energy of an impurity in a 
condensate is shown as a function of the coupling constant a, for different values of the temperature. 
From top to bottom, these are /? = 100,50,20,10,8,6,4,2, corresponding in our example to T = 
3.5, 7, 18, 35, 44, 58, 88 and 175 nK, respectively. The dashed curve shows the strong-coupling Pekar 
approximation. The inset zooms in on the small a region: the dashed line in inset shows the second 
order perturbation result. 



/? = 2, where thermal depletion of the condensate starts to be appreciable. Standard cooling 
schemes, as mentioned, can cool down to roughly /5 = 15. The dashed line shows the weak- 
coupling perturbative result valid for temperature zero (and without any wave number cut- 
off). As predicted by Timmermans and co-workers, at small a -C 1 and /3 ^ 1 the polaronic 
contribution to the energy is positive. It reaches a maximum around a = 1.5 — 2.0 and then 
decreases. The polaronic energy contribution becomes negative (indicative of a transition 
from an unbound state to a self-trapped state) at a critical coupling strength etc — 3 for 
T — > 0. This critical value goes to zero for increasing temperatures. The perturbative 
solution is seen to fit well at low coupling. At larger coupling, the dashed curve indicates 
the strong-coupling variational result (^^: it is significantly larger than the result of the 
full variation with M and Q free parameters, indicating that a Gaussian wave function may 
not be as suitable as it is for Frohlich polarons in the solid state. We emphasize that the 
polaronic energy calculated here is the contribution from ([6]) and does not contain the terms 
Egp + NqVib{0) which have a known dependence on the various tunable parameters Nq, ajB 
and asB, and which complicate the experimental determination of the polaronic energy 
contribution. To observe polaronic effects, it may be more straightforward to measure the 
shift in effective mass of the impurity. 
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FIG. 2: (Color online) The variational mass parameter M is shown as a function of a. The effective 
polaron mass is to a good approximation given by rupoi = mj{l + M/mi). The curves correspond 
to different temperatures {(3 = 100,50,20,10,8,6,4,2, from top to bottom), at cut-off Kc = 200. 
The inset shows the effect of cutoff {Kc = 100, 200,400, 1000, from top to bottom), at /3 = 10. 



C. Effective mass increase 



The effective polaron mass nipoi can be derived from the path-integral propagation of a 
particle from r(0) to a nearby point r(r) = r(0) + vt by casting the resulting transition 
amplitude in the form exp{— mpo^f ^r/2}. Feynman notes that this procedure gives a value 
for TUpoi which is always within a few percent of 1 + M with M the variationally optimal 
mass M of the trial model [H, 2^. Figure 2 shows the result for M as a function of a, 



for different values of the temperature. For small values of the coupling strength, the mass 
increases linearly with a as predicted by perturbation theory. However, near a ~ 3.5, the 
behavior changes and the mass increases rapidly. 

The low-energy Bogoliubov excitations are sound waves, with a dispersion qualitatively 
similar to acoustic phonons. The effective mass of acoustic polarons - electrons interacting 
with acoustic phonons - shows a jump of several orders of magnitude in the effective mass (cf. 
[l3| ) above a critical coupling strength. In the present case we also note a faster increase in 
the mass above a critical coupling strength, even though the interaction amplitude is different 
from that between acoustic phonons and electrons. In the dilute atomic gas, the transition 
is much less dramatic, and becomes smoother as temperature increases. We believe the 
smoothness of the crossover is not an artifact of the path-integral formalism, since in the 
case of the acoustic polarons it is the same formalism that predicts a discontinuous jump. For 
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acoustic polarons, it is also predicted that the sharpness of the transition depends strongly 
on the cutoff: above a critical value of Kc a discontinuity appears in the mass as a function 
of a. So it is worthwhile to study the dependence of M on the value of a cut-off to the k 
integrations in fll8l) also in the present case. We find that increasing the cutoff sharpens the 
transition (as can be seen in the inset of Fig. [2]), but no discontinuity arises as it does for the 
acoustic polaron. The difference is due to the fact that although the Bogoliubov excitation 
dispersion and the interaction amplitude show a similar /c-dependence as in the case of 
acoustic polarons, this is only true in the limit <C 1, and large deviations already occur 
for /c ^ 1, the relevant length scale of the problem. Yet even though the transition is not 
as abrupt as for acoustic polarons, it is possible to distinguish two regimes. In conjunction 
with the crossover from positive to negative values of the free energy, this is again indicative 
of a transition between an quasi-free (unbound) impurity and self-trapping for the impurity. 



D. Polaron radius 

The second variational parameter, f2, is linked to the polaron radius. Within the model 
system described by the action functional iSq, expression (fT3|) . the expectation value of the 
square of the relative coordinate for the impurity - boson mass system is given by 

(r ) = -7^ coth — 28) 

The square root of this is a measure of the localization length of the impurity wave function. 
For strong coupling, this expectation value converges to the expectation value with respect 
to the variational wave function formulated by Landau and Pekar[27|. 

Fig. ([3]) shows the polaron radius a/ (r^) as a function of the coupling strength, for 
different values of /3, at a cut-off Kc = 200. As the coupling strength is increased, the 
polaron radius decreases, indicating a stronger confinement of the impurity wave function. In 
particular the large a behavior might still depend on the cut-off, although a 1/a dependence 
is expected independently of Kc- Near a ~ 3 — 3.5, we find that a local maximum develops 
and the behavior of the polaron size as a function of a becomes non-monotonic. It is not 
clear whether this local maximum is an intrinsic feature of the BEC-impurity polaron, or 
whether it is an artifact from the particular variational model used. 



V. CONCLUSION 

When the Bogoliubov approximation applies, the Hamiltonian describing the impurity 
in a condensate can be cast in the form of the Frohlich polaron Hamiltonian. The physics 
becomes similar to that of a Frohlich polaron, where for the impurity in the BEC the Bogoli- 
ubov excitations have taken the role of the phonons and the interaction strength is related 
to the impurity-boson and boson-boson scattering lengths. The most accurate description of 
polaron physics in the case of electron-phonon interactions is given by Feynman's variational 
treatment, which moreover allows to study the temperature dependence of observables such 
as the effective mass and the free energy. We have applied a Jensen- Feynman path-integral 
type calculation to the case of the impurity in a condensate and derived expression f|T8l) for 
the free energy. Both in the free energy and the effective mass, a critical value of the coupling 
strength a ~ 3.5 can be identified where the system crosses over from the weak-coupling 
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FIG. 3: (Color online) The polaron radius in units of the healing length is plotted as a function of the 
coupling constant a, for different temperatures. From top to bottom, (3 = 100, 50, 20, 10, 8, 6, 4, 2. 
Around a ~ 3 nonmonotonous behavior develops as a function of a when the temperature is 
reduced. 



to the strong-coupling regime. These regimes show a qualitatively different behavior of the 
effective mass and free energy. The sharp increase in the the effective mass in the strong- 
coupling regime hints at a transition from an quasi-free state to a large-mass state similar to 
that for acoustic polarons. It has been pointed out 13| that a transition from the quasi-free 
regime to the large-mass regime is impossible to realize in most semiconductors and III-V 
compounds, even in alkali halides. However, the present results indicate that it might be 
attainable in ultracold gases. For this purpose, we investigated the experimentally relevant 
values of the system parameters, and derived a useful expression (1251) relating a to the var- 
ious parameters in the case of a condensate in the Thomas-Fermi regime. This opens up 
the prospect to reach and investigate the strong-coupling regime in ultracold gases, whereas 
this regime has hitherto be inaccessible in the solid state. 
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